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ABSTRACT 

We present an analysis of the X-ray spectra of two strongly magnetic cataclysmic 
variables, DP Leo and WW Hor, made using XMM-Newton. Both systems were in 
intermediate levels of accretion. Hard optically thin X-ray emission from the shocked 
accreting gas was detected from both systems, while a soft blackbody X-ray component 
from the heated surface was detected only in DP Leo. We suggest that the lack of a soft 
X-ray component in WW Hor is due to the fact that the accretion area is larger than 
in previous observations with a resulting lower temperature for the re-processed hard 
X-rays. Using a multi-temperature model of the post-shock flow, we estimate that the 
white dwarf in both systems has a mass greater than 1 Mq . The implications of this 
result are discussed. We demonstrate that the 'soft X-ray excess' observed in many 
magnetic cataclysmic variables can be partially attributed to using an inappropriate 
model for the hard X-ray emission. 

Key v^rords: Stars: binaries: eclipsing - Stars: magnetic field - Stars: novae, cata- 
clysmic variables - Stars: individual: DP Leo, WW Hor - X-rays: stars 
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1 INTRODUCTION 

Polars, or AM Her systems, are interacting binary systems 
in which the magnetic field of the accreting white dwarf 
is strong enough to prevent the formation of an accretion 
disc. It also forces the spin period of the white dwarf to be 
synchronised with the binary orbital period. In these sys- 
tems, the accreting material leaves the red dwarf secondary, 
follows a ballistic trajectory before getting coupled by the 
magnetic field of the white dwarf and funnelled directly onto 
the white dwarf surface. 

Just above the white dwarf photosphere, the accretion 
flow is strongly shocked, at which point the temperature 
reaches -^SOkeV. Some of this radiation is intercepted by 
the white dwarf where it is thermalised and re-radiated as 
soft X-rays. An additional source of soft X-rays is dense 
blobs of material which penetrate into the white dwarf pho- 
tosphere and heat it directly. In addition, there is a cyclotron 
component which is strongest in the optical band (for the 
magnetic field strength found in polars: ~10-200MG). This 
arises from electrons spiralling around the magnetic field 
lines close to the shock. 

Prior to the launch of the ROSAT X-ray satellite in 



1990 the relative strength of the three emission components 
was a major issue (eg King & Watson 1987). Observations 
of several polars showed that the strength of the soft X-ray 
component relative to the hard X-ray component was in ex- 
cess of that predicted by simple accretion models (eg Heise 
et al 1985). With the results derived using ROSAT it be- 
came clear that many polars showed this 'soft X-ray excess' 
(Ramsay et al. 1994) . During epochs when the accretion rate 
was in an 'intermediate' level of accretion, Ramsay, Cropper 
& Mason (1995) found that the relative strength of the soft 
X-ray flux was less than in 'high' accretion states: indeed 
AM Her showed a 'soft X-ray deficiency' in this state. They 
hypothesized that for lower accretion rates, a shock does not 
form over the whole of the accretion region. 



A significant uncertainty in the results of Ramsay et al 
(1994) was the poor constraint on the hard X-ray flux since 
ROSAT did not extend to energies much higher than ~2keV. 
With the launch of XMM-Newton we are able, for the flrst 
time, to detect both the soft and hard X-ray components 
clearly. This paper presents an analysis of the spectra of the 
first polars to be observed using XMM-Newton. 
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Object 



Observation 
Date 



Eflfective exposure (ksec) 
PN, MOS 



DP Leo 22 Nov 2000 

WW Hor 4 Dec 2000 



19, 18 
21, 23 



Table 1 . The observation log of the XMM-Newton observations 
of DP Leo and WW Hor. In the case of the DP Leo observations, 
time intervals of enhanced particle background were excluded. 



EPIC detector 



response file used 



PN 


epn_fs20.Y9.thin.rmf 


MOSl 


ml_thinlv9ql9t5r4_all_15.rsp 


MOS2 


m2_thinlv9ql9t5r4_all_15.rsp 



Table 2. The response files used in the analysis. 



2 OBSERVATIONS 

DP Leo and WW Hor were observed using XMM-Newton in 
Nov and Dec 2000 respectively: the observation log is shown 
in Table Py. Both objects were detected in all three EPIC 
detectors (Turner et al 2001, Stiider et al 2001), which were 
configured in full window mode and used the thin filter. 
Neither object was detected in the RGS detectors. During 
the DP Leo observations, there were time intervals of en- 
hanced particle background. These intervals were removed 
from the spectral analysis, reducing the exposure by 16ksec 
in the EPIC detectors. The observations of WW Hor were 
not significantly effected by increased particle background. 
By comparing the X-ray and optical fiux with previous ob- 
servations, it was concluded that both systems were in an 
intermediate state of accretion (cf Pandel et al 2001). 

Before extracting source spectra, the data were pro- 
cessed using the first public release of the XMM-Newton 
Science Analysis System. Spectra were extracted from the 
EPIC PN camera data using an aperture of ~ 19 in radius 
centered on the source, chosen so that the aperture did not 
cover more than one CCD. This encompasses ~85 percent 
of the integrated PSF (Aschenbach et al 2000). In the case 
of the EPIC MOS data, the source was in the centre of the 
CCD and a radius of ~ 24 was used (encompassing ~90 
percent of the integrated PSF). Background spectra were 
extracted from the same CCD on which the source was de- 
tected, scaled and subtracted from the source spectra. 

Since the response of the EPIC PN detector is not well 
calibrated at present below ~0.2keV, energies below this 
were ignored in the subsequent analysis. Also the PN re- 
sponse function used assumes only single pixel events, so we 
extracted only these events to construct our PN spectrum. 
For the EPIC MOS detectors, which are currently better cal- 
ibrated at lower energies, we fitted energies above O.lkeV. 
The response files used in the analysis are listed in Table El 



3 THE MODEL FOR THE X-RAY SPECTRUM 

The hard X-ray component produced in the hot post-shock 
flow has traditionally been modelled using a single temper- 
ature thermal bremsstrahlung component. We have devel- 
oped a physically more realistic model of the post-shock flow 
(Wu, Chanmugam & Shaviv 1994, Cropper, Ramsay & Wu 



1998), which takes account of its multi-temperature nature 
and the fact that a signiflcant fraction of the cooling takes 
place in the form of cyclotron radiation. It also includes the 
hard X-rays reflected from the surface of the white dwarf. 
In addition. Cropper et al (1999) included the effect of the 
varying gravitational force over the height of the post-shock 
flow. In this paper we use the model of Cropper et al (1999) 
to model the hard X-ray component and a blackbody to 
model any soft X-ray component. In addition, we account 
for absorption using a neutral absorption model. 



4 DP LEO 

The integrated EPIC PN and MOSl spectra of DP Leo are 
shown in Figure hi A soft X-ray component is seen together 
with a hard component - this is the first time that DP Leo 
has been detected above 2keV. Line emission from Fe Ka is 
detected at ~6.7keV. 

We fitted all three EPIC X-ray spectra separately us- 
ing the model described in 93. We fix the specific accre- 
tion rate, m, at 1 g s~^ cm" (appropriate for a system 
in an intermediate accretion state), the ratio of cyclotron 
to bremsstrahlung cooling, ts, at 10 (appropriate for the 
magnetic field strength seen in DP Leo: 31MG Cropper & 
Wickramasinghe 1993) and the metal abundance at solar. 
The fit is not sensitive to the exact values chosen for these 
parameters. 

The best fit parameters are shown in Table 0: good fits 
are achieved. The absorption towards DP Leo is low, < 10^" 
cm~^ and the temperature of the blackbody is similar to 
that found for other polars (eg Ramsay et al 1994). To de- 
termine the implied mass of the white dwarf we assume the 
Nauenberg (1972) mass-radius relationship for white dwarfs. 
This yields a white dwarf mass greater than I.SMq. 

We also show in Figure 1, the spectra centered on 7keV. 
Our model fits underestimate the line flux for the Fe Ka 
emission lines. To investigate this further, we varied the 
metal abundance, es and m. No solution was achievable that 
was consistent with both the line and continuum flux. We 
address this further in §H. 

We show in Table H^the unabsorbed, bolometric fluxes 
for the soft and hard X-ray components. The soft X-ray 
luminosity is deflned as Lsoft,boi = '^fsoft.hoiaec{i — P)cP , 
where we assume that the soft X-ray emission is optically 
thick and can be approximated by a small thin slab of mate- 
rial, the unabsorbed bolometric blackbody flux is fsoft,boi, d 
is the distance, i is the system inclination and /3 is the angle 
between the accretion region and the spin axis. For DP Leo, 
i = 80° and /3 = 100° (Bailey et al 1993), while Biermann et 
al (1985) find d > 380 pc. Because we have been able to ex- 
tend to lower energies in the MOS detectors compared to the 
PN detector (since they are currently better calibrated) , the 
absorption is more constrained using the MOS detectors. In 
determining the flux from the soft component using the PN 
data we have constrained the absorption to be < 9 x lO'^'' 
cm~'^ - the upper limit derived from the MOS detectors. 

We define the hard X-ray luminosity as Lhard,boi = 
'^Ttfhard,boid^, where fhard,boi is the unabsorbed, bolometric 
flux from the hard component. Since a fraction of the hard 
X-ray flux is directed towards the white dwarf and some of 
that is reflected back towards the observer, we switch the 
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reflected component to zero after the flnal fit to determine 
the intrinsic flux from the opticaUy thin shocked material. 

We also show in Table p^ our estimate of the ra- 
tio Lsoft,boi/Lhard,boi- Using data obtained using ROSAT, 
Ramsay et al (1994) were able to put an upper limit only 
on this ratio since the hard X-ray component was not sig- 
niflcantly detected. Now using the XMM-Newton data, we 
can place relatively good constraints on the ratio because we 
sample both the hard and soft X-ray components. If there 
were no dense blobs of material in the accretion stream, 
we expect Lsoft,boi/ Lhard,boi ~ 0.5 (in some polars in a 
high state of accretion, it is found that this ratio is much 
greater than this; on the other hand, for systems in interme- 
diate states of accretion, the ratio can be much lower). The 
Lsoft.boi/Lhard,boi ratio in DP Leo, is around unity, imply- 
ing consistency with the standard Lamb and Masters (1979) 
model. 

Cropper, Wu & Ramsay (2000) suggest that using a sin- 
gle temperature thermal bremsstrahlung to model the hard 
X-ray component will contribute to the apparent 'soft X- 
ray excess' seen in many polars. This is because a fraction 
of the photons attributed to the soft X-ray component actu- 
ally originates in the lowest levels of the post-shock flow 
rather than from re-processing by the white dwarf. Our 
multi-temperature shock model intrinsically emits a signif- 
icant amount of soft X-rays at the base of the post-shock 
flow. It is interesting to determine the soft to hard X-ray 
ratio if we assume an absorbed blackbody plus single tem- 
perature thermal bremsstrahlung model as used by Ramsay 
et al (1994). Fixing the thermal bremsstrahlung tempera- 
ture at 30keV (as there), and adding a Gaussian to model 
the Fe Ka line emission at ~6.7keV, we find that the ratio, 
Lsoft,boi/Lhard.boi is twicc the value determined using the 
multi-temperature model. This suggests that when compar- 
ing the relative strength of the shocked and re-processed 
X-ray components, it is essential to use the most physically 
realistic model of the post-shock flow available. 



5 WW HOR 

The integrated EPIC PN and MOSl spectra of WW Hor 
are shown in Figure ti. Whilst we observe hard X-ray emis- 
sion above 2keV in WW Hor for the first time, there is no 
evidence for a soft blackbody component with the sort of 
temperature seen in DP Leo. This is unusual, since even for 
systems in a reduced rate of accretion a soft X-ray compo- 
nent is expected. Weak Fe Ka line emission is present at 6.7 
keV. In the PN spectrum there is evidence for an emission 
line due to Fe fluorescence near 6.4keV. 

We modelled the integrated EPIC spectra from each de- 
tector separately using the model described in §H together 
with a neutral absorber. As in DP Leo we flx rh, at 1 g s~^ 
cm~'^, es at 10 and the metal abundance at solar. We ob- 
tain good flts: the best flt model parameters are shown in 
Table 0. We flnd a very low absorption column. The mass 
of the white dwarf is ~1.1Mq. Using a distance of 430pc 
(Bailey et al 1988) we estimate the unabsorbed, hard bolo- 
metric luminosity to be Luard = 1 — 2 x 10^^ ergs s~^. This 
is similar to that found for other AM Her systems in inter- 
mediate accretion states (Ramsay, Cropper & Mason 1995) . 
Indeed, when we compare the phase averaged observed flux 



in the 0.1-2.0keV band in XMM-Newton (2.1 x 10"^^ ergs 
s"^ cm-^in MOSl) with the ROSAT observation of WW 
Hor, made when it was in an intermediate state (2.3 x 10~^^ 
ergs s~^ cm^^), we flnd that they are very similar. 

We can investigate whether a low temperature soft X- 
ray component could be absorbed by a modest absorption 
column. In the lowest energy bands the best calibrated in- 
strument is the EPIC MOSl detector. We added a black- 
body component with temperature kTi,i,—10 and 15 eV to 
the model used above (cf kTbb ~3Q eV for DP Leo). We find 
that such a low temperature blackbody is 'hidden' by an ab- 
sorption column consistent with our upper limit (5.6x10^^ 
cm~'^). The corresponding values of L^oft is 8.5x10^° and 
4.3x10^" ergs s~^for 10 and 15eV respectively. This corre- 
sponds to Li,oft/Lhard=0.4: aud 0.2: these values are broadly 
similar to that found for polars in intermediate accretion 
states (Ramsay, Cropper & Mason 1995). 



6 DISCUSSION 

When WW Hor was observed using ROSAT it was found 
that both a blackbody and a bremsstrahlung component 
were required to fit the X-ray spectrum (Ramsay, Cropper 
& Mason 1995). When WW Hor was observed using XMM- 
Newton, it was at very similar X-ray flux levels compared 
to the ROSAT observation. Further, the X-ray bolometric 
luminosity of WW Hor is approximately similar to DP Leo 
(implying similar mass transfer rates) if we assume similar 
distances. Why do we therefore observe a soft X-ray compo- 
nent in DP Leo and not WW Hor? 

We suggest that the accretion region was more spread 
out in the XMM-Newton observation of WW Hor compared 
to the ROSAT observation, thereby decreasing the local 
heating rate and causing the temperature of any soft X- 
ray component to decrease. If it was sufficiently low (say 
fcT=10 or 15eV), we showed in §ra that even a very modest 
amount of absorption would prevent us from detecting this 
soft X-ray component. It is not clear why the accretion re- 
gion would cover a larger area of the white dwarf. One pos- 
sible reason could be that for relatively low mass transfer 
rates, accretion occurs onto a wide range of magnetic az- 
imuths, resulting in a more extended accretion region. This 
is thought to be one of the reasons why intermediate polars 
(which have less strong magnetic fields than polars) do not 
show a soft X-ray component (eg King & Lasota 1990). 

We now briefly comment on the derived masses for the 
white dwarfs in these systems. In the case of WW Hor, we 
find a mass of --l.lA'/eand for DP Leo >1.3M0(table |). 
The mass of the white dwarf in both systems is at the upper 
end of expectations. There are no other accurate estimates 
for the mass of the white dwarf in these two systems. Several 
magnetic CVs were found to have similar masses to WW Hor 
in the survey of such systems using RXTE data and the same 
emission model as used here (Ramsay 2000). Similarly, the 
mass of Sirius B is 1.034±0.026 (Holberg et al 1998). Other 
massive non-accreting white dwarfs have been found (eg RE 
J0317-858; 1.31-1.37Mq, Ferrario et al 1997). In the case of 
DP Leo, however, the mass of the white dwarf is very heavy 
indeed and would be the heaviest white dwarf within an ac- 
creting magnetic system yet found. Ramsay (2000) found 
that the mass distribution of isolated non-magnetic white 
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channel energy (keV) 

Figure 1. The EPIC PN (upper), MOSl (lower) X-ray spectra of DP Leo. The soUd hnes show the best fit using an absorbed blackbody 
plus multi-temperature accretion shock model. 




channel energy (keV) 



Figure 2. As for figure 1, but for WW Hor 
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(xlOis (eV) 
cm~2) 



Ml 
(Mq) 



FluXao/t,boi 

(xl0~^2 erg 
s^-"- cm~2) 



Lsoft.bol 

(xlO^l 
ergs s"-'-) 



(xlO"'^^ erg 



^hard,bol 

(xlO^i 
ergs s"-') 



Lsoft.boll 
^hard,bol 



PN 

MOSl 

MOS2 



0.0 

1.4" 



-1.4 



1.04^ 



1+4.2 
-1.2 

0.0+1-6 



1+0.05 
-0.06 
-I lo+O.lO 
J^-J^^-0.09 

-,+0.06 



0.99" 



■0.05 



7 6+°-'* 
'•"-0.3 

+ 1.0 
1.5 

+0.4 
0.5 



12.8 
7.6 



1.7 

2 

1.7 



+0.1 
0.1 
+0.2 
-0.2 
+0.1 
0.1 



(dof) 



DP Leo 


PN 

MOSl 
MOS2 


0.0+20 

9 ,+6.9 
^•1-0.6 

q q + 5.7 

•^•■^-1.8 


23+t 
30-6 
28I5 


1.40-0.12 
1.40-0.09 
1.40-0.08 


(5 y+36.2 
"•'-3.7 

^•5+^j 


3.4+l»,6 
1.3+^-'' 

1 fi+2-O 


44+0.8 
^•^-0.5 
c c+1.0 

^•6_i.o 

c c+1.3 

^•6Io.9 


1 0+0.1 

^■•^-0.2 

1 0+0.2 
^■•^-o.i 


49+29 

^•"-0.5 

1 2+2-'' 
-^•^-0.8 


1.12 (126) 
0.88 (90) 
1.05 (72) 



1.25 (123) 
1.02 (116) 
0.98 (165) 



Table 3. The best fit parameters to the EPIC PN, MOSl and MOS2 data of DP Leo and WW Hor using the multi-temperature model 
described in 33. A distance of 400pc is assumed for DP Leo (Biermann et al 1985 find d > 380pc) and 430pc for WW Hor (Bailey et al 
1998). In determining the fluxes and luminosities of DP Leo using the PN data, we have chosen an upper limit of Nu; < 9 X lO^^ cm~2 
which is the upper limit derived from the MOS data (which has a better response at lower energies). 



dwarfs and white dwarfs in magnetic CVs are significantly 
different: the white dwarfs in magnetic CVs are biased to- 
wards heavier masses. However, it was not clear if this was 
due to selection effects: it is possible that high mass polars 
are more likely to be discovered than low mass polars. 

In calculating the mass of the white dwarf we use 
the Nauenberg (1972) approximation for the white dwarf 
mass-radius relationship: this approximates the Hamada & 
Salpeter (1961) mass-radius relationship for carbon white 
dwarfs over the range 0.3-1.2 Mq, and is consistent with 
the best mass-radius data for the white dwarfs 40 EriB 
and Sirius B (Provencal et al 1988). However the appro- 
priate mass-radius relationship is still a matter of debate. 
In the case of DP Leo, where we estimate a white dwarf 
of mass >1.3M0, there is some additional uncertainty for 
our mass estimate since for the very heaviest white dwarfs 
the mass-radius relationship is very uncertain. Hamada & 
Salpeter (1961) also show the mass-radius relationship for Fe 
core white dwarfs. White dwarfs with Fe cores are thought 
to be formed when a white dwarf with a ONeMg core ig- 
nites through an accretion-drived collapse (Isern, Canal & 
Labay 1991). The maximum mass for such a white dwarf 
is ~1.1M0. Using this mass-radius relationship we obtain a 
mass of ~1.1M0, again close to the maximum. 

Alternatively, it is possible that additional cooling pro- 
cesses not included in our model could become important for 
white dwarfs greater than I.OMq. This would have the ef- 
fect of modifying the temperature-height profile of the post- 
shock region and hence the implied X-ray spectrum. If the 
volume of gas emitting ^lOkeV was underestimated, then 
this may also account for the deficit of Fe line emission in DP 
Leo (§M). In general, however, additional cooling terms work 
to soften the X-ray spectrum, resulting in a higher gravi- 
tational potential (and hence white dwarf mass) to fit the 
observed spectrum. It is also possible that the absorption is 
more complex than a simple neutral absorber, as was found 
to be the case in the polar BY Cam (Done & Magdziarz 
1998). This would make the continuum appear harder, so 
overestimating the shock temperature and hence the white 
dwarf mass. The higher temperature continuum would also 
underpredict the plasma iron lines fluxes. Given that both 
these features are seen in our model fits, we added a partial 
covering component and also an ionized absorption compo- 
nent to our model. The resulting fits and masses were not 
significantly improved with the exception of the MOSl data 
which gave masses which were 0.1-0.2 MQlighter. Data of 



higher signal to noise would be required to test the necessity 
of more complex absorption models. 
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